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We analyze antiferromagnetism and superconductivity in novel Fe— based superconductors within 
the weak-coupling, itinerant model of electron and hole pockets near (0, 0) and (n, tv) in the folded 
Brillouin zone. We discuss the interaction Hamiltonian, the nesting, the RG flow of the couplings 
at energies above and below the Fermi energy, and the interplay between SDW magnetism, su- 
perconductivity and charge orbital order. We argue that SDW antiferromagnetism wins at zero 
doping but looses to superconductivity upon doping. We show that the most likely symmetry of 
the superconducting gap is Aig in the folded zone. This gap has no nodes on the Fermi surface but 
changes sign between hole and electron pockets. We also argue that at weak coupling, this pairing 
predominantly comes not from spin fluctuation exchange but from a direct pair hopping between 
hole and electron pockets. 
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I. INTRODUCTION 

Recent discovery of superconductivity (SC) in the iron- 
based layered pnictides with Tc ranging between 26 and 
52K generated enormous interest in the physics of these 
materials, which also hold hold strong potential for ap- 
plications. — . SC has been observed in oxygen containing 
1111 systems RFeAsO, where R=La, Nd, Sm, Pr, Gd, in 
oxygen-frcc 122 systems AFC2AS2, where A=Ba, Sr, Ca, 
and in several other classes of materials like LiFeAs with 
111 structure and a-FeSe with 11 structure —i^i^i^iii^. 

In several respects, the pnictides are similar to the 
cuprates. Like the cuprates, the pnictides are highly two- 
dimensional, their parent materials show antiferromag- 
netic long-range order below 150 K^i^i^°'^^i^^ , and super- 
conductivity occurs upon doping of either electrons^ii^i^ 
or holes^ into FeAs layers. This lead to early con- 
jectures that the physics of the pnictides is similar to 
that of the cuprates and involves insulating Mott behav- 
jQj.^i3ji4ji5ji6 j^esistivity measurements, however, showed 
that iron pnictides remain itinerant down to zero dop- 
ing, although the jury is still deliberating whether some 
signatures of Mott physics have been observed at higher 
energies (see e.g., -'-^1-'-''' ). Other evidences for itinerant 
behavior include 

• a relatively small value of the observed magnetic 
moment per Fe atom in the magnetically ordered 
phase - 12 - 16% of 2/ii3 in nil materialsi^ii^. 

• a good agreement between electronic band struc- 
ture calculationsi^iiai20,2i^^ ^^^^ ARPES and 

magneto-oscillation measurements of the Fermi 
surface (FS) and electronic states 

• Drude-likc behavior of the optical conductivity at 
small frequencies^ 

Although there are certain variations in the crystal 
structure between different classes of pnictides, the low- 



energy electronic structure is likely the same for all sys- 
tems and consists of two small hole pockets at the center 
of the Brillouin zone (BZ) and two small electron pock- 
ets centered around M points (Fig[T]). The M points are 
located qi = (0,7r/a) and q2 = {n/a,0) in the unfolded 
BZ (one Fe atom in the unit cell) and at identical points 
ki ~ (7r/a, 7r/a) and k2 = {n/a, —n/a) in the folded BZ 
(two Fe atoms in the unit cell, a = a-\/2) (see Figl2|). 
The relations between the momenta in the folded and un- 
folded zones are = {qx + <}y)/V2, ky = {q^ - qy)/^/2. 
The unfolded BZ includes only Fe states, the folded BZ 
(the correct zone) takes into account the fact that only 
a half of Fe states are actually identical in pnictides be- 
cause of As which resides either above or below Fe plane 
(Figl^l). Below we will use the folded BZ. Throughout 
the paper we define M point as Q = (7r/a,7r/a) and set 
a = 1. 

The goal of this paper is to summarize recent work 
done in collaboration with D. Efremov and I. Eremini^ 
and by myself on the weak-coupling, Fcrmi-liquid analy- 
sis of magnetic and superconducting instabilities in the 
pnictides. I will address several issues: 

• what interactions cause magnetism and SC? 

• are SC and magnetism competing orders? 

• is charge order possible? 

• what is the symmetry of the SC gap and how to 
overcome intra-pocket repulsion 

• is the pairing magnetically mediated? 

Central to weak coupling analysis is the idea of a near- 
nesting between electron and hole pockets at zero dop- 
ing. The nesting does not mean that any of the FS 
has parallel pieces but rather implies that, by moving a 
hole FS by Q, one obtains a near-perfect match with an 
electron FS. This FS geometry is in a reasonable agree- 
ment with ARPES and magneto-oscillation measure- 
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FIG. 1: (color online) A simplified geometry of Fe— based su- 
perconductors used in the present work. The Fermi surface 
consists of an electron pocket around (tt, n) (black solid cir- 
cle), and a hole pocket of roughly equal size around (0, 0) (blue 
solid circle) . A near-perfect nesting between hole and electron 
pockets means that, by moving a hole FS by {n,n), one ob- 
tains a near-perfect match with an electron FS. Upon electron 
doping, the size of the electron pocket increases (dashed blue 
black), what breaks the nesting. Upon hole doping, the 
size of a hole FS increases, what again breaks the nesting. 
-|-A and —A are the values of the superconducting gaps on 
the two FS for an s+ superconducting state. (From Ref.— .) 
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FIG. 2; (color online) Folded and unfolded BZ (courtesy of I. 
Eremin). Unfolded BZ is a x a square (one Fe atom per unit 
cell). Unfolded zone (the correct one) is constructed to take 
into account the fact that only a half of Fe states are actually 
identical because As, whose projection onto an Fe plane is 
at a center of a square, is located slightly above or below an 
Fe plane. The folded zone has two Fe atoms in the unit cell 
and is a square with dimensions y/2a x \/2a, rotated by 45'^ 
compared to the unit cell in unfolded zone. 



ments 24, 25, 26, 27, 28, 29 _ j^g^g been known from the stud- 
ies of chromium and its alloys in the 70t h^^i'^^ and from 
generic theoretical studies of "excitonic insulators"— 
that such nesting leads to a spin-density-wave (SDW) 
order with momentum Q already at weak coupling in 
the same way as SC order appears at weak coupling in 



a BCS superconductivity. For pnictidcs, the role of nest- 
ing for SDW order was emphasized by Cvetkovic and 
Tesanovio^ and Barzykin and Gorkov^^. Upon doping, 
either by holes or by electrons, one of the two pockets 
gets relatively larger and the nesting breaks down. In 
this situation, SDW order gets weakeri^MLi^^i^^ and at 
some doping superconductivity emerges. How pnictides 
evolve from a magnet to a superconductor is not clear 
at the moment, and there are experimental evidences for 
both first order and continuous transitions. Theoreti- 
cally, both situations arc possible.— 

What interaction causes superconductivity and what 
is the symmetry of the superconducting gap are the two 
most intriguing issues for the pnictidcs. A conventional 
phonon-mediated s-wave superconductivity is an unlikely 
possibility because electron-phonon coupling calculated 
from first principles is quite small^. An electronic mech- 
anism is therefore more likely. Mazin et al^ conjectured, 
by analogy with the cuprates, that SC in pnictidcs is me- 
diated by antiferromagnetic spin fluctuations. This pair- 
ing mechanism is most effective if the superconducting 
gap changes the sign under the momentum shift by Q. 
For the cuprate FS, this unambiguously leads to dx2-y2 
superconductivity. For pnictides, the requirement is that 
the gap must change sign between hole and electron FS. 
This does not unambiguously determines the gap struc- 
ture, but most natural would be the gap which is a con- 
stant A along a hole FS and a constant —A along an 
electron FS. In the classification of the cigcnfunctions of 
the tetragonal D^h group, this gap belongs to Aig repre- 
sentation and is roughly A(k) = A(cosfca; + cos ky)/ 2^. 
We will refer to this gap symmetry as s^. The s"*" gap 
has been found as the most likely candidate in some of 
the RPA studies based on a 5-orbital Hubbard model^, 
a 2-band spin-fluctuation model^, and in the renormal- 
ization group (RG) analysis^i^ (see below). The s+ gap 
structure also emerges in the analysis based on localized 
spin models^. 

There are. however, two potential problems with the 
spin- fluctuation mechanism. First, the close proximity 
to a magnetic phase does not a'priori guarantee that 
the pairing is magnetically mediated. This is particu- 
larly true for pnictides because of two separate FS in 
which case there are multiple interactions between hole 
and electron states, and the interaction which gives rise 
to magnetism is not necessary the same interaction that 
gives rise to SC. This has to be verified in the calcu- 
lations. Wc will argue below that the full pairing in- 
teraction docs have a component which represents the 
exchange by soft dynamic magnetic fluctuations peaked 
at Q. However, we will argue that this component is 
subleading at weak coupling, and the dominant pairing 
component is a direct pair-hopping term. This does not 
change the outcome, though, that the pairing is in the 
s'^ channel. 

Second, intra-pocket repulsion does not cancel out 
from the s'^ pairing problem because average gap along 
either hole or electron FS is non-zero. The intra- 
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pocket interaction term in the Hamiltonian is quite likely 
stronger than the pair-hopping ter m^^i^^ such that at 
this level 3+ channel is repulsive. One has to see to 
where the interactions flow at small energies and whether 
the renormalized pair-hopping term eventually becomes 
larger than intra-pocket repulsion. 

This flow of couplings is outside RPA in which the in- 
teractions have the same values as in the Hamiltonian.— 
Whether or not s+ superconductivity is favored in RPA 
then likely depends on the values of the bare interac- 
tions in the underlying 5-orbital model (intra-and inter- 
orbital Hubbard interactions, intra-orbital exchange and 
pair-hopping terms). While some researchers found the 
s+ stated, others^ found that more likely candidates 
are two nearly degenerate states in which the gap has 
nodes on one of the FS sheets and no nodes on the other. 
One such state is another extended s-wave state with 
A(k) « Acos^cos^, the second is d^^^y^ state with 

A(k) sa A sin ^ sin In the unfolded BZ, these two are 
cosq^+cosqy and cosq^^—cosqy, respectively (see Refi^). 
These two states transform into each other under the 
translation by Q, but are decoupled in the gap equation. 
The extended s-wave state with A(k) oc cos ^ cos ^ is 
not orthogonal to A(k) cx cosk^ + cosfcj, as both are 
members of the Aig representation of D/^h group. How- 
ever, for small hole and electron pockets, the matrix ele- 
ment between cos + cos ky and cos kx/2 cos ky/2 states 
is small, and the two can be treated "almost" indepen- 
dently. 

The gaps with A(fc) cx cos cos ^ and A(fc) cx 

sin -Y sin -j- arc less sensitive to intra-pocket repulsion 
as the average of A(fc) along one of the FS sheets is 
zero, and win within RPA over cosk^ -\- cosky state when 
intra-pocket repulsion is the strongest interaction. On 
the other hand, these two states are less favorable candi- 
dates once we include inter-pocket attraction because the 
magnitude of one of the gaps at the FS is small, of order 
kp. Other pairing states, like d^y states in the folded BZ 
with A(k) = A sin fca: sin fcj, have also been proposed^ 
and arc favored in a parameter range in which the pair- 
ing interaction is peaked at momenta smaller than 2k p. 

From experimental perspective, the issue of the gap 
symmetry is also unsettled. ARPES^^i2^>^i21 and An- 
dreev spectroscopy^ measurements have been inter- 
preted as evidence for a nodeless gap, either a pure s- 
wave gap or an s"*" gap. The resonance observed in neu- 
tron measurements below T^r^ is consistent with an s+ 
gg^p39^ but not with a pure s— wave gap or cos ^ cos ^ 

or sin ^ sin ^ gaps. On the other hand, nuclear mag- 
netic resonance (NMR) data— i^i^ and some of the pen- 
etration depth dat a^'^1^^ were interpreted as evidence for 
the nodes in the gap. Some, but not all of these data can 
still be reasonably fitted by a model of an s+ SC with or- 
dinary impurities^i^i^'^'^. Several group o^^'^^i^° sug- 
gested specific measurements which could potentially un- 
ambiguously distinguish between different pairing sym- 
metries. 



In the bulk of the paper we discuss weak-coupling ap- 
proach to the pnictides. The most essential results of our 
calculations are 

• magnetism and SC emerge due to the interplay be- 
tween pair hopping and inter-pocket forward scat- 
tering 

• magnetism and SC are competing orders; mag- 
netism wins for perfect nesting but looses to SC 
upon doping 

• orbital charge order appears at T only slightly be- 
low magnetic ordering temperature 

• the pairing symmetry is s"*", intra-pocket repulsion 
is not an obstacle, at least when kp is small 

• pairing interaction is not magnetically mediated - 
it predominantly comes from a direct pair hopping 
between hole and electron FS 

Weak-coupling approach is certainly only a first step 
in the understanding of the physics of the pnictides and 
should be followed by more involved calculations which 
include the dynamics of the interactions. These dynamic 
calculations should allow one to understand what sets 
the upper cutoffs for the effective interactions in density- 
wave and pairing channels, to obtain the numbers for 
T/v and Tc, and to study quantitatively the behavior of 
various obscrvablcs. Only after that one should be able 
to settle the issue whether weak/moderate coupling ap- 
proach to the pnictides captures most of the physics, or 
strong coupling effects associated with Mott physics are 
also essential. 

The paper is organized as follows. In Sec|TT] we in- 
troduce the model and discuss the approximations. In 
Scc lIIII we discuss renormalization group flow first at en- 
ergies above Ep and then at energies below Ep. We show 
that RG equations are different in the two regions. In 
Sec lIVI we discuss density- wave and pairing instabilities 
and address the issues whether magnetism and super- 
conductivity are competing orders, can a charge density- 
wave order emerge, and is there a symmetry between 
different orders. In Sec|V] we discuss whether the pair- 
ing can be viewed as mediated by spin fluctuations. In 
Scc lVIl we present the conclusions. 



II. THE LOW-ENERGY MODEL 

We model iron pnictides by an itinerant electron sys- 
tem with two hole pockets centered at (0, 0) and two elec- 
tron pockets centered at Q = (tt, tt) in the folded BZ. We 
assume that the interactions do not distinguish between 
the two hole FS (and between the two electron FS) and 
restrict with one hole and one electron FS adding combi- 
natoric factors to the interaction effects whenever neces- 
sary. There are fine magnetic effects associated with two 
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rather and one hole (and electron) FS^i^ but we will 
not discuss them here. 

There are three key assumptions of our weak-coupling 
approach to the pnictides. 

• We assume that the interactions are small com- 
pared to the bandwidth such that the low-energy 
physics is determined solely by fermions with mo- 
menta near (0,0) and Q. The alternative descrip- 
tion in terms of lattice spin models (like Ji — J2 
modeli^iii^iii'^'^) is valid if the interaction is com- 
parable or stronger than the bandwidth. 

• We assume that there is near-nesting between elec- 
tron and hole pockets (by moving a hole FS by Q 
one obtains a near-perfect match with an electron 
FS). Nesting does not have to be perfect, but a typ- 
ical energy associated with non- nesting should be 
smaller than the Fermi energy Ep- 

• We assume that both hole and electron pockets are 
small compared to the size of the BZ. In energy 
units, this implies that the Fermi energy is much 
smaller than W. 

The last two assumptions are at least partly consistent 
with ARPES and magneto-oscillation experiments and 
with band theory which all show that there is a near- 
nesting between at least one hole and one electron pock- 
ets, and also that both pockets are small and the Fermi 
energies (top of the band for hole pocket and bottom of 
the band for electron pocket) are of order O.leV, much 
smaller than the bandwidth W, which is roughly 2eV. 
The smallness of Ep compared to W is important for 



our description because it opens up a relatively large en- 
ergy window between W and Ep in which, on one hand, a 
low-energy description of nested electron and hole pock- 
ets is still valid, and, on the other hand, it is unimportant 
where precisely a fermion resides near one of the FS. We 
show below that in this particular regime, particle-hole 
and particle-particle channels of fermionic interaction be- 
come degenerate, and all five relevant couplings flow log- 
arithmically under parquet RG. 

The actual situation is indeed more complex than 
our simlified model. In particular, ARPES data show 
that there are two roughly equivalent electron pock- 
ets, but only one small hole pocket of about the same 
size^'^i^'^. The other hole pocket is larger in size, 
and the superconducting gap along this pocket is about 
two times smaller than the gap along the smaller hole 
pocket. Theoretical calculationsi^iiSi2ai^i2^i^i^ also 
show a somewhat more complex electronic structure, 
with another, fifth band in a close proximity to a Fermi 
level. 

We label fermions near (0, 0) as c— fermions and 
fermions near Q as /-fermions. We assume that they 
have a parabolic dispersion 
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For convenience, we shift momenta of /—fermions by Q 
such that = — Cp. In 2D, fermionic density of states 
is energy independent and we will use this in the RG 
analysis. 



There arc five different interactions involving c- 
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The momentum conservation is assumed in all terms. 

We label the couplings with subindcx "0" to emphasize 
that these are bare couplings. The terms with vf'' and 
are intraband interactions, the terms with uf'"' and 
^72*^^ are interband interactions with momentum transfer 

and Q, respectively, and the term with is inter- 
band pair hopping. We show these interactions graph- 
ically in FigO We approximate all five interactions by 
momentum-independent constants because of small sizes 
of hole and electron FS. 

Below we will be using dimensionless interactions 

uT^uf^N, (3) 
where Aq = m/{2T:) is the fermionic density of states. 



Weak-coupling approach implies that u\ <1. 

Eq. ^ is the effective Hamiltonian for the interaction 
between low-energy fermions. It can be obtained, in prin- 
ciple, from the underlying five-orbital model^ with intra- 
orbital and inter-orbital Hubbard interactions U and V 
and the Hund's rule coupling J. The relations between 
our and U, V, J should, however, be rather involved 
because all five Fe orbitals contribute to low-energy elec- 
tronic states. The actual situation in pnictides is even 
more complex because Fe states also hybridize with As 
p— states^. 

Note that all vertices in Eq. ^ are (5-functions in spin 
indices, i.e., all bare interactions are in the charge chan- 
nel. Bare spin-spin interaction terms with spin matrices 
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FIG. 3: Five difFerent interactions between low-energy 
fermions. Solid and dashed lines represent fermions from c— 
band (near k = 0) and f — band (near k = Q = (tt, tt)). (From 
Ref.^.) 



in the vertices are also possible in principle, but at weak 
coupling such terms are generally smaller than density- 
density interactions. Spin vertices, however, appear once 
the interactions get dressed up by singular particle-hole 
bubbles (see below). 



III. RENORMALIZATION GROUP ANALYSIS 

Interactions between low-energy fermions may gener- 
ally give rise to density-wave and pairing instabilities of 
a Fermi liquid. These instabilities can be handled within 
weak-coupling theories only if the corresponding response 
functions are strongly enhanced, and the enhancements 
overcome the smallness of the interactions. One example 
of such behavior is a pairing instability for which one only 
needs an attraction in the corresponding pairing channel 
- a well known fact is that an arbitrary small attrac- 
tion already makes the system unstable against pairing. 
Mathematically, this is due to the fact that at T = 0, 
particle-particle polarization bubble Ilpp(q,n) made out 
of two c— fermions or out of two /—fermions diverges 
logarithmically when q,fl^O: 



(27r)3 



iVolog 



where G''^^^ {k,uj) = {i 



A 



max {n, vpq) 



(4) 



^ frec-fermion prop- 
agator and A is the upper cutoff. 

Density-wave instabilities generally require Ui to be 
above a threshold of order one and are not captured 
within weak-coupling theory. Pnictides are special in this 
regard because of the nesting. The nesting implies that a 
particle-hole polarization bubble np^(g-|-Q, f2) with mo- 
mentum transfer near Q, made out of one c— fcrmion and 



one /— fermion, also diverges logarithmically at g, J7 — > 0: 
lil{{q + Ct,n)^ [ ^^G^,{k,u;)Gl{k + q + Q,iU + n) 
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(5) 



Particle-hole bubbles determine the responses of a sys- 
tem to density-wave perturbations in spin and charge 
channels, and the divergence of n^{(Q,0) imphes that 
the system becomes unstable towards a particular 
density-wave order already at a weak coupling, once the 
interaction in the corresponding channel becomes attrac- 
tive. If there is more than one channel with attractive 
interaction, the order first appears in a channel where 
the attraction is the strongest. 

The obvious issue then is what are the values of the 
interactions in different channels. In a generic weakly 
coupled Fermi liquid, the interactions are the combina- 
tions of the parameters of the Hamiltonian. Pnictides, 
however, are again special because, as we said, Ep is 
much smaller than the bandwidth W. We will see that 
a conventional weak-coupling behavior, in which differ- 
ent channels compete, only holds a.t E < Ep. At larger 
energies Ep < E < W the system flows towards a novel 
fixed point at which all five interactions tend to infinity, 
but their ratios tend to universal values. We will see be- 
low that at this point, the symmetry extends to 5*0(6) 
(Refi^) , and the pairing channel, the SDW channel, and 
also the orbital charge density-wave (CDW) channel be- 
come completely equivalent and critical. 

Such symmetry, however, is only present right at the 
fixed point, where all five renormalized couplings are in- 
finite, i.e., the system is on the verge of an instability. 
The issue then is whether the fixed point behavior is ac- 
tually reached already above Ep. This depends on initial 
conditions. 

If the fixed point is reached at some E > Ep , then the 
instability of a Fermi liquid is determined by an SO (6) 
fixed point, and the spin and the charge density- wave 
orders and the pairing order emerge simultaneously, as 
the components of a six-dimensional order parameter. In 
this situation normal state Fermi liquid behavior does 
not exist at energies smaller than Ep. Furthermore, de- 
viations from a perfect nesting should not change the 
system behavior as long as such deviations do not affect 
energies above Ep. If the fixed point is not reached at 
E > Ep, then the system flow at Ep < E < W sets 
the values of the renormalized, but still finite couplings 
at E ~ Ep. These renormalized couplings then act as 
the "bare" couplings for the theory at energies E < Ep, 
at which spin density- wave, charge density- wave and su- 
perconductivity become competing orders, each develops 
at its own energy below Ep. Deviations from a perfect 
nesting are obviously more relevant in this situation as 
smaller energies are involved. 

In the pnictides, both magnetic T/v and superconduct- 
ing Tc are substantially smaller than Ep ^ lOOOA', so 
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very likely density-wave and pairing instabilities emerge (a) 

from energies below Ep, where SO {6) symmetry is bro- Pi P3 Pi ^ ^ P3 

ken. Below we assume that this is the case and consider K _ \ \ 

separately the flow of the couplings above Ep and below /} ~ ) ) 



E 



F- 



A. Energies larger than Ef 

In a generic weak-coupling case wp"* < 1, integrating 
out fermions with energies between Ep and W does not 
substantially affect the system because the corrections to 
the couplings contain higher powers of u\^\ However, in 
a situation when Ep << W such renormalizations con- 
tain uf^^ \ogW/ E, where i? is a running scale between 
Ep and W. As a result, interactions flow logarithmically 
from their bare values at energies of order W to their 
renormalized values at a scale E. When uf'' log W/E 
become of order one, the renormalized couplings may dif- 
fer substantially from the parameters of the Hamiltonian. 
Without nesting, such logarithmic renormalization would 
only occur in the particle-particle channel. When nesting 
is present, the renormalizations in both particle-particle 
and particle hole channels are logarithmic. 

We emphasize that, at energies larger than Ep, the 
renormalization of the vertices from Eq. ^ is indepen- 
dent on what their total and transferred momenta are 
compared to kp, e.g., a particle-particle vertex with any 
total momentum k < kp of two c— or two /—fermions 
undergoes the same renormalization. Similarly, vertices 
with transferred momenta between c— and /—fermions 
k' = Q + q undergo the same renormalization for all 
q kp. 

As it is customary to weak coupling theories, below we 
only consider renormalizations which contain powers of 
uf''^ log W/E and neglect regular renormalizations which 
bring additional powers of u^^\ In a diagrammatic lan- 
guage, this amounts to summing up series of diagrams 
and keeping the highest power of the logarithm at every 
order. The presence of two orthogonal channels of loga- 
rithmic renormalizations implies that one needs to sum 
up parquet series of diagrams (see, e.g. Refi^). We ex- 
plicitly verified, by evaluating diagrams up to third order 
in uf'\ that the system is renormalizable, i.e., that, in- 
stead of summing up infinite series of diagrams, one can 
write up one-loop parquet RG equations. The deriva- 
tion of the RG equations is straightforward (see Fig. [4]). 
Evaluating second-order diagrams and collecting combi- 
natoric pre- factors, we obtain the set of equations 

ill = uf + 

U2 = 2U2(U1 - U2) 

U3 = U3(4ui - 2u2 - U4 - U5) 

2 2 

U4 = — U3 — M4 

"5 = -W3 - "s- (6) 
where the derivatives are with respect to L = 




FIG. 4; Diagrams for the one-loop vertex renormalizations. 
The diagrams for the renormalizations of iti and U3 terms are 
shown. The diagrams for the renormalization of U2, 114, and 
U5 terms are obtained in a similar fashion. (From Ref.— .) 

(1/2) logW/E. The factor 1/2 reflects the fact that each 
of the fcrmionic bands extends to energies above Ep only 
in one direction. This factor was neglected in Ref.—. 
Similar, though not identical equations have been ob- 
tained in the weak-coupling analysis of the cuprates for 
"t— only" dispersionSft. 

The RG equations for the couplings U4 and M5 are 
identical and preserve particle-hole symmetry. Below we 
only use U4. Note that all couplings remain momentum- 
independent under one-loop RG. The momentum depen- 
dence appears beyond the leading logarithmic approxi- 
mation, and is small at weak coupling. 

We see from Eq. ([6]) that the pair hopping term U3 is not 

generated by other interactions, i.e., U3 = if Ug"' = 0. 
When uf^ = 0, the set ^ decouples into 

ill = 

(2u2 - ill) = -(2u2 - wi)^ 

iiA = -ul (7) 
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FIG. 5: (color online) The RG flow of Eqn. ([6]) in variables 
U4,/u3 and u\/uz (from Ref.— ). The flxed point is U4/U3 — 
— 1/\/5, u\/u-j, = \l\fh. For definiteness, we used boundary 
conditions m^°'' = u\ \ itg"' = O.Im^'^'', and set Mq = 0. (From 
Ref.^.) 



In this special case, particle-hole and particle-particle 
channels decouple: M4 = uf^ j (\ + ^.hu^^^ \ogW j E) is 
renormalized in the particle-hole channel and flows to 
zero for uf' > 0. [By the same reason, strong Coulomb 
repulsion renormalizes down at low energies in conven- 
tional SC and allows electron-phonon interaction to over- 
come it at small frequencies^]. The renormalization 
of the inter-pocket interaction ui comes from particle- 
hole channel, is of opposite sign, and ui — u^^^ /{I — 
O.^u^i"^ logW/E) increases when E decreases. The renor- 
malization of 2m2 — wi also comes from particle-hole chan- 
nel, and U2 tends to ui/2 when 2^2° "* — u'f' > 0. 

Once u^^^ is finite, the system moves into the basin of 
attraction of another fixed point at which all couplings si- 
multaneously tend to infinity and U3 becomes the largest. 
Indeed, introducing U2 = aui, U3 = and U4 = 7U1, 
we can re-write Eq. (|6]) as 

a — uia (1 — 2a — , 

/3 = ui/3(3-2(a + 7)-/32), 

7 = -ui(l+7) (/3'+7)- 

(8) 

The new fixed point corresponds to d = /3 = 7 = 0, i.e., 
to 

/3(3-2(a-f 7)-/3') =0, 
(1 + 7) + 7) = 0. 

(9) 

For /? 7^ (i.e., U3 7^ 0), the solution of the set is a = 
0,7 = -1, /3 = ±V5. There are no other solutions for 



/? 7^ 0. Expanding near this fixed point we obtain 
ill ~ 6u^, a — — 4ui a, 

X bx ui [2a + 2z + x), z = — 4iti z. (10) 

where 2 = 1+7 a; = — 5. Solving this set, we 
find that ui increases, while a, x and z all scale as "^^^ 
and vanish when ui diverges. This implies that the fixed 
point specified by ([9]) is a stable fixed point. 

We see that \u^\ becomes the largest coupling at the 
fixed point, while U4 = — wi, and ui — \u^\/\/^. The sign 
of W3 at the fixed point is the same as the sign of Ug"^ , i.e., 
it is positive for u"^^ > and negative for Ug"' < 0. Using 

a (X we find that near a fixed point U2 oc jugl"'^/'^. 

We also verified that a fixed point with U3 ~ (i.e., 
/? = 0) is an unstable fixed point, while the one for which 
U3 is the largest coupling is a stable fixed point. 

In Figl5]we show the RG flow obtained by a numerical 
solution of Eq. ^ . We see that the system indeed flows 

towards a stable flxed point for any Ug"'' 7^ 0. On the 
other hand, this stable fixed point is reached only at an 
infinite coupling. As long as Ui are finite, the system 
flows towards the stable fixed point, but does not reach 
it. A scale when Ui diverge is log W/E < 1/u*, where 
u* is some linear combination of the bare u^. Like we 
said, we assume that at E = Ep-, the couplings are still 
flnite, i.e., this fixed point is not yet reached. In practical 
terms, this implies that logW/Ep < 1/u*. 

How the system fiows depends on the values and signs 
of the bare Ui. We assume that all bare interactions are 
repulsive, Ui > 0. Then, under RG flow, ui and U3 in- 
crease, while M4 decreases, passes through zero, changes 
sign, becomes negative and approaches —ui. The cou- 
pling U2 also increases, but becomes relatively small com- 
pared to other couplings. The sign change of 114 is the 
consequence of the coupling between particle-hole and 
particle-particle channels. If particle-hole channel was 
not logarithmic, the RG equation for ug would reduce to 
Mg = — 2ugU4. Solving this equation together with the 

equation for U4 one would obtain that, for uf^ > |ug'^''|, 
the intra-pocket repulsion U4 would preserve its sign and 
just decrease logarithmically under RG together with the 
pair-hopping term. 

B. Energies smaller than Ef 

The RG flow described by Eqs ^ stops atE ^ Ep- At 
smaller energies the independence of the RG equations of 
the total and transferred momenta is lost, and, e.g., the 
vertex ug with zero total momentum of two c— fermions 
is renormalized differently from the same vertex with the 
total momentum of order kp. To put it simply, only two 
types of vertices continue to flow logarithmically at ener- 
gies below Ep - the vertices with zero total momentum 
of two c— or two /—fermions, and (for the case of a per- 
fect nesting) the vertices with the momentum transfer 
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between c and /— fcrmions exactly equal to Q. The ver- 
tices with zero total momentum are M4 and M3 terms, the 
vertices with the momentum transfer Q are ui, U2 ;and 
M3 terms. The vertex U3 with zero total momentum gen- 
erally has an arbitrary momentum transfer Q + (7, where 
q ^ kp- The vertex U3 with the momentum transfer Q 
has the total incoming momentum of order kp. For all 
other vertices, the RG flow is cut at Ep. 

We will use the values of Ui at Ep as the "bare" cou- 
plings for the theory at £^ < Ep and label them as 
{ui = Ui{E = Ep)). The upper limit for the RG at 
E < Ep is then obviously Ep. We label the couplings 
with zero total momentum as 1x3(0) and ^4(0) and the 
couplings with momentum transfer Q as ui{Q), U2{Q), 
and U3a(Q), and u^biQ)- The separation into u^a and 
M3h is due to the fact that there arc two different W3 
vertices with momentum transfer Q: the coupling u^a 
corresponds to ^3 = pi + Q, in the U3 vertex in Fig. [31 
while U3b corresponds to pi = ^4 -I- Q. 

The RG equations for the relevant couplings are ob- 
tained in the same way as at higher energies, but now 



one has to carefully analyze which diagrams still yield 
log Ep/E and in which diagrams the logarithm is cut. 
For example, the renormalization of 1x3(0) comes from the 
first diagram in Fig. the renormalization of u^aiQ) 
comes from the second, third, and fourth diagrams, and 
the renormalization of ustiQ) comes only from the third 
diagram. Evaluating the diagrams, we obtain 



dui{Q) _ 2 



dusbiQ) 



dL 
du2{Q) 
dL 

dU3a{Q) 

dL 
dusjO) 
dL 



uiiQ)+ul,{Q), -^'"^ ^2u3biQ)uiiQ), 

= 2u2{Q)iui{Q) - U2iQ)) + 2u3a{Q){u3biQ) - U3aiQ)), 
= -4u3a(Q)M2(Q)) + 2u2iQ)u3biQ) + 2u3aiQ)ui{Q), 



_2«3(0)u4(0), ^^^^ul{0)~ul{0). 



(11) 



where now L = log Ep/E. From (jlip we immediately 
obtain 



^ (uiiQ) + u3biQ)) - iu,{Q) + u3biQ))\ ^ iu3b{Q) - u,{Q)) 



iu3b{Q)-u,iQ))^, 



-JY (uiiQ) + U3b{Q) - 2{U2{Q) + U3a{Q))) = {ui{Q) + U3b{Q) - 2{u2{Q) + U3a(Q)))', 

dL 

4t iMQ) - U3b{Q) - 2{U2{Q) - U3a{Q))) = {UliQ) - U3b{Q) - 2{u2iQ) - «3a(Q)))', 

dL 



^ (W3(0) - H4(0)) = {U3{0) - U4(0))2, ^ (U3(0) + U4(0)) 



-(M3(0) + M4(0))2 



These equations can be easily solved and yield 

Ml -I- U3 



UliQ) + usbiQ) 



1 - (Ml + U3) log ^ ' 



UliQ) + U3biQ) - 2(u2(g) + U3a{Q)) = 

UliQ) - "36(g) - 2iu2iQ) - U3aiQ)) 
U3 - Ui 



M3(0)-U4(0) = 



1 - (U3 - "4) log ^ 



UliQ) - u3biQ) = 

1 - iui 

Ui-U3~ 2U2 

1 - (iii - U3 - 2u2) log ^ ' 

Ui + U3 — 2U2 

l-iui+U3-2u2)\or^''' 



Ui - U3 



U3 



)log^' 



U3(0) + U4(0) = 

I 



= E 
U3 + Ui 



1 + (U3 + Ui) log ^ ' 



(12) 



(13) 



We see that there are six different energies E at which 
different combinations of the couplings diverge or may 
diverge 



IV. DENSITY- WAVE AND PAIRING 
INSTABILITIES 

We now relate the divergences of these six combina- 
tions of the couplings with density-wave and pairing in- 



stabilities for our model. We searched for SDW and 
CDW instabilities with momentum Q and with cither 
real or imaginary order parameter, and for a SC insta- 
bility cither in pure s channel (the gaps Ac and A / have 
the same sign), or in s+ channel (the gaps Ac and A/ 
have opposite sign). The instabilities with momentum- 
dependent order parameter, like a nematic instability^ 
or a d— wave superconductivit y"^^'^^ do not occur in our 
model because we set bare interactions to be momentum- 
independent and considered only the leading logarith- 
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+ 



-<r >U4(0) + 



M - 




^SDW^y^^ ^SDW 






FIG. 6: Diagrammatic representation of the renormalized ver- 
tices in SDW, CDW, and SC channels. The shaded triangles 
are hilly renormalized vertices, double lines represent running 
couplings given by Eq. p3p . 



mic renormalizations of the vertices. The momentum 
dependence of the vertices comes from subleading, non- 
logarithmic renormaUzations which we neglected. 

Density-wave and pairing susceptibilities are obtained 
by standard means: by introducing an infinitesimal cou- 
pling in a particular channel and evaluating the response. 
For this, we add to the Hamiltonian three extra terms 

k 
k 

fe k 

(14) 

with complex Agdw, ^cdw, and real A^/, and evaluate 
how these terms are renormalized by Ui . The correspond- 
ing diagrams are presented in Fig. [6l In analytic form, 
the fully renormalized vertices at an energy E (or, equiv- 
alently, at a temperature T) are given by 



A{- = A,(i + r,iog§: 



(15) 



where j corresponds to either SC, or SDW, or CDW or- 
der parameters. We label corresponding Tj as F^'^^^y, 

^CDW^ and t'qq \ where r, i mean real or imaginary 
density-wave order, and s, s"^ mean s— wave or extended 
wave SC order, respectively An instability towards a 
particular SDW, CDW or SC order occurs at a tempera- 
ture Tg^yy, "^CDW^ '^^ -^sc ' which the corresponding 
Fj(i?) diverges. 



FIG. 7: (color online) The values of effective interactions Fj 
in various density-wave and superconducting channels at an 
energy E ~ Ef vs L — log VK/-Bf. The effective interac- 
tions are linear combinations of Ui (Eq. (HHJ), which are the 
solutions of the RG set (|6]). For definiteness we set bare pa- 



rameters Mj"' 



,(0) 



Mo, and 



m(0) 



O.lMo. AH 



Fj are in units of uq. Observe that the strongest interaction 

is in SDW channel, and that f"sc changes sign and becomes 
attractive once L exceeds some critical value. (From Ref.— 
with permission from the authors. ) 



From Fig. [SI we have 



SDW 
CDW 
CDW 



UliQ) 

:Ul(Q) 

-UliQ) 



U3biQ), 



r 



SDW 



2(U2(Q) 
2(W2(Q) 



UliQ) - usbiQ), 

U3a{Q)), 
USaiQ)), 



^sl = -("4(0) + U3(0)), r^sc' = "3(0) - W4(0) (16) 

We see that these vertices contain exactly the same 
six combinations of parameters which we analyzed in the 
previous section. Combining Eqs. ([T^, and (fTBl) we 
find that all six channels are decoupled at energies below 
Ep, all Tj satisfy the same RG equation 



dF, 
dL 



= n 



(17) 



The interplay between different instabilities is then de- 
termined by the bare values of Tj in different channels, 
which are Tj{E ^ Ep)- This imphes that SDW, CDW, 
and SC orders are competing orders in this approxima- 
tion. The first instability occurs in the channel for which 
the coupling aX E ^ Ep is the largest. 

At i? ~ Ep, Ui = Ui, and we have {Tj{E ~ Ep) = Tj) 



(r) 

SDW 



(r) 

CDW 



Ui 

- Ul 



"3, 



f(i) 
^ SDW 

2u2, T 



Ul 



^SC = "4 



CDW 

"3, r^",^' =U4-U3 

We see that Tgj^^y diverges when 

1 



Ep 
^og — 



Us, 

= Ul 



1 



U3 



SDW "3 + Ul 



(18) 
(19) 
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if U3 + ui > 0. The vertex T^^^^ diverges at 
Ef I 1 



log' 



E r 



SDW 



if ui > U3. The vertex diverges 



at 



1 

log — 



-pi 

^ CDW 



(21) 



if wi + Us — 2u2 > 0. The vertex ^cdw diverges at 
Ef I 1 



log- 



E 



CDW Ui-U2-2U2 

if Ml — -03 — 2u2 > 0. The vertex diverges when 
Ef 1 -1 



(22) 



log 



E r 



if M3 + U4 < 0. The vertex Tip diverges when 



log 



Ep 
~E 



1 



1 



r 



sc 



U3 - Ui 



(23) 



(24) 



if 'O3 > M4. 

Strictly speaking, only the largest E (for which 
log Ef/E is the smallest) has physical meaning within 
this approach as a Fermi liquid does not exist below this 
energy. However, when two or more critical E are close, 
the system may have a sequence of phase transitions. 

In Fig. [7] we show how these Tj evolve depending on 
the value oi L = logW/EF- We see that for our choice 
of M-" > ^sFiw t'^^ largest, and the highest critical E 
is when F^^;^ diverges. This implies that for a perfect 
nesting the first instability of a normal state occurs in the 
SDW channel and leads to a real SDW order parameter. 
This happens at 



SFiw « exp 



Ml + M3 



(25) 



At the same time, the energies at which Tcd\y and F|^ 
diverge are only slightly smaller if the bare couplings Ui 
are close to their fixed point values U3 = v^^i, M4 = 
—Ml, M2 << Mi: instability towards CDW order with an 
imaginary order parameter (an orbital order) occurs at 



Thow cxp 



1 



Ml 



M3 — 2m2 



and the instability in the s+ SC channel occurs at 
Tsc « exp ' 



M3 - M4 



(26) 



(27) 



Whether CDW and SC orders emerge as additional or- 
ders at a smaller T requires a separate analysis as density- 
wave and the pairing susceptibilities change in the pres- 
ence of an SDW order. Ref. — found that additional 
orders do not occur for a perfect nesting. 



We emphasize again the the special role of the RG 
renormalization at energies higher than Ef- The bare 



(20) couplings Mg"^ and M4"'' are the parameters of the Hamilto- 
nian, and without RG flow at intermediate energies, there 
would be no divergence of the couplings at zero total mo- 



,(0) 



mentum if u 



(0) 

'(0) 



> M 



(0) 



Moreover, it is very likely that 



in pnictides m^'^'' is larger than Mg*^^ because 
from the Hubbard repulsion within a given orbital, while 
Mg*^^ comes from inter-orbital processes which is generally 
though to be weake r^^i^^i^^ . However, we know that the 
renormalization from energies above Ef strongly affects 
M4 and not only reduces its magnitude but forces it to 
change sign and become negative. In this situation, at 



,(0) 



comes 



(0) 



,(0) 



< 0, 



energies of order Ef, M3 — M4 > even if u 
and the fully renormalized M3(0) — M4(0) diverges at E 
given by (pi)) . 

When nesting becomes non-perfect (e.g. upon dop- 
ing), the interplay between different instabilities changes. 
The pairing channel still remains logarithmic, and Tg^ 
is given by ([?7|) . In the density- wave channels, the po- 
larization bubble made of c— and /— fermions with mo- 
mentum transfer Q remains logarithmic {log Ef/ E) only 
for E larger than the energy scale 5 associated with non- 
nesting. At smaller energies, the logarithm is cut. Ac- 
cordingly, TgJJ^r{S) and TQjjyy{5) decrease, and at some 
deviation from a perfect nesting, Tgij^r{5) = T§q. At 
larger S, the first instability of a Fermi-liquid is into a SC 
state with an s"*" symmetry of a superconducting gap. 

The actual transformation from an SDW to a SC state 
is more involved because SDW state for a non-perfect 
doping eventually becomes incommensurate^^*^, like in 
chromium^2iS. The incommensurate state is a magnetic 
analog of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
stated of a superconductor in a magnetic field^i^. De- 
pending on whether or not a SDW order becomes incom- 
mensurate before superconductivity emerges, the transi- 
tion from a SDW antiferromagnet to a 5+ SC is either 
first order or involves an intermediate phase in which 
SC order co-exists with an incommensurate SDW order— 
(see Fig. [5]). A peak in the spin susceptibility at an in- 
commensurate momentum near Q. has been observed 
in the weak-coupling, RPA analysis of a multi-orbital 
Hubbard- like model^i (five Fe bands). This and other 
studies^i^ also reported smaller peaks in the suscepti- 
bility at a smaller momenta, possibly related to 2kF- 



A. SO{6) symmetry 

For completeness, we also consider a situation when 
hole and electron pockets are so small in size that the 
regime W > E > Ef extends down to an energy at which 
Ui diverge, and the ratios of the couplings approach their 
fixed point values M3/M1 = -\/5, U4 = — mi and M2/M1 
0. In this situation, the system becomes unstable even 
before it enters a regime E < Ef- One can immediately 



see from (1161) that in this situation F' 



SDW 



CDW 
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5 / 23ir s / 27ir 



FIG. 8: (color online) Phase diagram for Ts/Tc = 1.5 (a) 

and Ts/Tc = 5 (b), where T^ = rj^„ and Tc = T^c- S 
measures a deviation from a perfect nesting in energy units 
(from Ref.— ). On the left panel thick solid and dashed lines 
represent second-order SDW and SC transitions, dotted line 
represents first order transition between commensurate SDW 
(SDWo) and s+ SC states, and light lines denoted by 5a and 
Sm are instability lines of SC and SDW phases, surrounding 
a first-order transition line. On the right panel, an incom- 
mensurate SDW order appears below Ts{5) once it becomes 
smaller than T* = 0.56Ts > Tc- Below Tc, a new mixed phase 
appears in which incommensurate SDWg order co-exists with 
SC. At small T, there is no SDW, state without superconduc- 
tivity. The transitions into the mixed state are second order 
from a SC state and first order from a commensurate SDW 
state. The transition from the normal state (N) to SDW, 
state is second order (dashed-dotted line) , and from SDW, to 
SDWo is first order (dotted line). (From Ref.— with permis- 
sion from the authors.) 

, i.e., SDW order, orbital CDW order, and s+ SC 
order appear simultaneously. 

The system behavior at such critical point has been 
discussed in detail by Podolsky, Kee, and Kim^. They 
demonstrated that the fixed point Hamiltonian has 
50(6) symmetry, and found 15 operators which form the 
Lie algebra of the 50(6) group. They found that the 
three interaction channels are indeed equivalent, and the 
order parameter is a six-component vector. Two com- 
ponents are the pairing amplitude and the phase, three 
components are spin projections, and the remaining com- 
ponent is an imaginary charge order parameter. The di- 
rection of this six-component vector is not specified by 
the Hamiltonian and set up by by a spontaneous break- 
ing of the 50(6) symmetry. An even larger symmetry is 
expected once one takes into consideration the fact that 
there are two hole and two electron FS^ 

In pnictides, the temperatures of magnetic and super- 
conducting instabilities are 10-20 times smaller than Ep, 
and it is therefore unlikely that the normal state becomes 
unstable already at £^ > Ep, although to firmly establish 
this would require more involved calculations than the 
weak-coupling analysis which we present here. At the 
same time, even approximate 50(6) symmetry should 



afi'ect system properties at intermediate energies. For 
example, a dispersing magnetic resonance mode in an s+ 
superconducting state has a velocity v = u_f/\/2 (Refj^), 
the same as the velocity of the Anderson-Bogolubov 
phase mode in a 2D uncharged superconductor. A num- 
ber of collective properties associated with an approxi- 
mate 50(6) symmetry were considered in—. 



V. IS THE PAIRING MAGNETICALLY 
MEDIATED? 

In this section we consider in more detail the case when 
an instability of a normal state occurs well below Ep. We 
found earlier that the flow towards 50(6) fixed point is 
cut at Ep, and at smaller energies SC, SDW, and CDW 
channels decouple, i.e., SC, SDW, and CDW become 
competing orders. Each develops independently, and the 
one with the highest instability temperature wins. 

This is true however only in the leading logarithmic ap- 
proximation. Beyond it, there is still a residual coupling 
between SC and density-wave channels, and the issue is 
whether this residual coupling may lead to a new physics. 
In particular, near a SDW instability, magnetic suscepti- 
bility near Q is strongly enhanced. If there is a residual 
coupling between SC and SDW channels, the pairing in- 
teraction should generally have a component which can 
be viewed as the exchange of antiferromagnetic spin fluc- 
tuations. At weak coupling, this component has an extra 
smallness simply because it comes from residual inter- 
action, but this smallness may be compensated by the 
enhancement of the antiferromagnetic spin susceptibility 
Xs(Q)- If the enhancement of Xs(Q) wins over the overall 
smallness of the residual interaction, spin fluctuation ex- 
change becomes the dominant component of the pairing 
interaction. 

The idea that the pairing near a magnetic instability 
is mediated by spin fluctuations has been applied to var- 
ious systems, most notably to the cupratcs^, where this 
mechanism of electronic pairing unambiguously yields 
d.j.2_y2 symmetry of the pairing gap. In the cuprates, 
however, there is only one large hole FS in the normal 
(not pseudogap) state, and the couplings are by no means 
small. 

Below we first derive the effective interaction and show 
that it indeed contains an extra term with the spin- 
fluctuation exchange. We then estimate the relative 
strength of this term compared to a direct 1*3 — M4 pairing 
interaction which we obtained earlier. 

For uf^ « 1, the residual coupling between SDW 
and SC channels can be actually obtained in a controlled 
way. The term which contributes to both SC and SDW 
is the pair-hopping term. Earlier in the paper we sep- 
arated this term into 1*3(0), which is the interaction with 
zero total momentum: 1*3(0) = u^ik, —k;p + Q, —p — Q), 
and u^aiQ) and u^b^Q) which are two interactions with 
momentum transfer Q: usa = u-i{k,p',k + Q,p — Q), and 
u^b = U3(k,p,p — Q,k + Q). The first two momenta 
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are for c— fermions, the other two are for /— fermions. 
Like before, we shift the momenta of /—fermions by 
Q. In these new notations, 1x3(0) — u^(k,—k;p,—p), 
U3a = U3{k,p]k,p), U3b = U3{k,p,p,k). Interactions u^a 
and u^b are density- wave vertices, and ^3(0) is a pairing 
vertex. Since both SC and SDW instabilities are confined 
to the FS, we restrict the interactions to |k| = |p| = kp- 

We now observe that, ah-cady at this stage, there 
is an overlap of measure zero between the pairing and 
density-wave vertices because the vertices u^a and u^t 
are defined for any angle between k and p, including 
p ~ — k. At this particular point u^a = u^^k, —k, ; k, — fc) 
and usb = uz{k,—k;—k,k) become components of the 
particle-particle vertex, albeit for a single momentum 
transfer. One should indeed exercise care to avoid dou- 
ble counting of the pairing interaction. In practice, this 
means that one has to subtract from u^a and u^b their 
bare value U3. Still, even after this subtraction, u^a and 
u^b diverge at a SDW transition, i.e., there exists a for- 
mally divergent component of the pairing potential, al- 
though at this stage it is only present for a single mo- 
mentum transfer. 

The overlap of measure zero has no physical effect, but 
we can extend the analysis by including into the pair- 
ing vertex the terms u^a and u^b at zero total momen- 
tum and at small but finite momentum and frequency 
transfers q, and oj. In our notations, these vertices are 
U3a{k, —k; k + q,—k — q) and u^bik, —k, —k — q,k + q) 
{q = (q, w)). These two vertices do not diverge at a 
SDW transition, but still are enhanced at small q. An 
important observation here is that these usa and u^b still 
come from the same series of diagrams as at g = - 
the only adjustment is in the argument of the logarithm. 
All other additions to the pairing interaction are regu- 
lar functions of q and can be safely neglected because of 
small overall factor. In other words, a potentially rele- 
vant pairing component due to residual interaction be- 
tween SDW and SC channels can be explicitly obtained 
from our earlier diagrammatic analysis, by extending it 
to a small but finite q and keeping the same diagrams. 
This reasoning parallels the one used in the derivation of 
the fully renormalized vertex functions in a Fermi liquid 
with predominantly forward scattering^. 

To obtain U3a{k, —k; k + q, —k — q) and U3b{k, — fc, — fc — 
q,k + q) we note that near a SDW instability, ui{Q) + 
usbiQ) diverge as (ui+us) / {1- {ui+U3)L) , while 
usiQ), U2{Q)~U3a{Q) and Ui{Q) + UsaiQ) - 2{u2{Q) + 
u^biQ) remain finite. Elementary calculations then show 
that 



notations 



U3a{k,-k,k 



-k-q) ^uJJiq) 



U3b{k, ~k, -k-q,k + q) = u^JJi^q) = ^ 



1 {ui+U3fL„ 



Al- {ui+U3)Lq 



1 [Ul + U3fLq 



2 1 - (Ui -I- U3)Lq 

(29) 

where Lq ~ logEp/Eq and Eq = E'^+avptl'^+bLJ^ , a,b ~ 
0(1). These expressions are only valid at small q, when 
I^f'^ and Ug^'^ are enhanced. Note the difference with a 



'■3a 



Coopcron in a disordered metal: there, Eq contains \uj\ 
term2^. 

Combining the two interactions given by (|29p with the 
original bare interaction in the pairing channel U3 — 
M4 we obtain irreducible, antisymmetrized pairing vertex 
Ta/S, 75(^1 —k,p, —p) in the form 

'^ap.,^s[k, -k,p, -p) = {U3 - UA)[5al35'^5 ~ SasSp^) 
+"3^''' {k - p)5ap5fS - u'JJ (fc - _p)(Sa5(S/37 (30) 

Using 6as5p^ = {SafsS^s + CTq^CT75)/2 and U3J{k - p) = 
Q.bu^J {k — p), we further obtain from (|30p 



^a[3,j5ik,-k,p,-p) = -{U3 ~ U4){5a[3SyS - ^aP^fS) 



ull^ (k ~ p)aat3ff^s (31) 



We see that the "original" U3 — U4 term has both charge 
and spin components, but the extra term has only the 
spin component and obviously can be interpreted as an 
exchange of dynamic antiferromagnetic spin fluctuations. 
Note the importance of the distinction between u^{/ and 



{u'^{/ = O.Sug^'^). If we didn't split U3 into two 
different components, the additional pairing term would 
have both spin and charge components. 

The term u'^l'^ [k — p) contains an extra power of u 
compared to the direct U3 — U4 term, but, on the other 
hand, is enhanced at small k—p. To see whether u'^l^ {k — 
p) is relevant for s"*" pairing, we need to average (j3ip over 
the FS to get s+ harmonic of the interaction F''^ =< 
F >i?5. Using |k— p| = 2k p smO/2 and integrating over 6 
we find that very near SDW transition angular integral is 
confined to small 9, where Eq. pQ]) is valid (c.f. RefiS). 
We then obtain right at the SDW transition 

ra^,75(fc, -k^p, -p) = -^(^3 - Ui){5al35^5 - CTapCT^s) 



-Cexp 



(Ul -I- U3) 



Ef 



(32) 



U3a{Q) ~ ^W36((3) 



1 {Ui + U3) 

4 1 - (wi + U3)L 



(28) 



Extending the result to a finite q, and subtracting the 
bare values to avoid double counting, we obtain, in our 



where C = 0{1). We see that the spin-fluctuation con- 
tribution to the pairing vertex is exponentially small at 
weak coupling, but contains 1/uj. To compare relative 
strength of regular and spin-fluctuation contributions, we 
take Lj ^ T and use for T the solution without spin- 
fluctuation contribution T = TJ^ ~ Epe-'^l^^^-^*'^. We 
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then obtain that at such T, spin fluctuation contribution 
contains extra 

exp — ( ^ — 7^ ^1 (33) 

At small coupling, this is an exponentially small factor. 
We see therefore that the spin-fluctuation contribution 
to the pairing is present, but at weak coupling it is ex- 
ponentially small compared to the regular term in the 
irreducible vertex, even if we bring the system to a SDW 
transition point. This implies that the residual coupling 
between SDW and SC channels does not change the sys- 
tem behavior at weak coupling — SDW and supercon- 
ductivity remain competing orders, each develops inde- 
pendent on the other. In other words, s+ pairing at weak 
coupling is not mediated by spin fluctuations but comes 
from a direct pair hopping from one Fermi pocket to the 
other. Spin fluctuation contribution well may become 
dominant at moderate/strong coupling, but this regime 
is beyond the scope of the present paper. 

Note the crucial role of the dynamics in our considera- 
tion. If the magnetically-mediate interaction was purely 
static, the spin-fluctuation term in (j32p would scale as 

U3b{0) and eventually become dominant near an SDW 
transition, when the divergence of u^b{0) overshadows 
the smallness of the overall factor. The presence of the 
dynamical component in magnetically-mediated pairing 
term keeps this term finite even when U3b(0) diverges. 

VI. CONCLUSIONS 

To conclude, in this paper we presented weak-coupling, 
Fermi liquid analysis of density- wave and superconduct- 
ing instabilities in Fe— pnictides. We modeled pnictidcs 
by a low-energy model of interacting electrons near small 
hole and electron FS located near (0, 0) and Q = (tt, tt) 
in the folded BZ. The interaction Hamiltonian contains 
five terms: intra-pocket repulsions U4 and U5 (which we 
set equal), inter-pocket density-density interaction terms 
Ml and U2, and a term U3 which describes a pair hop- 
ping from one pocket to the other. We assumed that the 
bare interactions uf''^ are positive and do not depend on 
a position of a fermion on either of the FS. 

We first considered energies between the bandwidth 
W ~ 2eV and the Fermi energy Ep ~ O.leV. We demon- 
strated that particle-hole and particle-particle channels 
are indistinguishable at these energies, and that all five 
couplings remain constants and logarithmically flow to- 
wards a fixed point with extended symmetry (identified 
as SO (6) in Refi^). In one- loop approximation which we 
used, the fiow is described by parquet RG equations. In 
the process of the RG flow, the pair-hopping term and the 
direct inter-pocket interaction term increase, while intra- 
pocket interaction M4 decreases, passes through zero and 
becomes negative below some energy scale. The interac- 
tion U2 also increases but becomes progressively smaller 
than other interactions. We argued that the RG flow 



favors superconductivity: even if the bare couplings are 
such that the bare pairing interaction is repulsive in all 
channels, the interaction in the extended s— wave pairing 
channel (s"*") becomes attractive once the renormalized 
M3 — U4 becomes positive. This definitely happens under 
RG because U4 eventually becomes negative. The order 
parameter in the s"^ channel is a constant on both hole 
and electron FS, but changes sign between them. 

The S0{6) fixed point is reached at an energy (temper- 
ature) when u* logW/E = 1, where u* is a linear combi- 
nation of five bare couplings. At this energy, all couplings 
diverge, and a Fermi liquid normal state becomes unsta- 
ble against the development of a six-component order 
parameter, whose (equivalent) components are SDW, or- 
bital CDW, and s+ SC. It would be extremely interesting 
if pnictides displayed this behavior, but this scenario is 
realized only if u* log W/ E = I occurs at E > Ep- 

Both magnetic and superconducting instabilities in the 
pnictides occur at T << Ep, and we argued that more 
likely scenario for this materials is that u* log W/ E < 1 
down to ~ Ep, and the instabilities are the results 
of further renormalizations at energies smaller than Ep. 
We considered the flow of the couplings at these ener- 
gies and found that the system behavior changes qual- 
itatively: SDW, CDW, and SC channels decouple, and 
each develops its own order at an energy (temperature) 
determined by the corresponding couplings at E = Ep. 
We found that for a near-perfect nesting, the instability 
at the highest energy is towards a conventional SDW or- 
der. When doping increases and the nesting becomes less 
perfect, the first instability eventually becomes a SC in- 
stability in the s"^ channel. The transition towards an or- 
bital CDW order (a CDW order with an imaginary order 
parameter) occurs only slightly below a SDW transition 
and, in principle, CDW fluctuations should be observ- 
able. 

Finally, we addressed the issue whether the pairing 
near a SDW instability can be viewed as mediated by spin 
fluctuations. To first order in the couplings, the pairing 
interaction comes directly from the pair hopping and has 
both charge and spin component. We explicitly demon- 
strated that the residual coupling between particle-hole 
and particle-particle channels at energies below Ep gives 
rise to an additional pairing component, which near an 
SDW instability reduces to the exchange by spin fiuctu- 
ations. However, we found that, even at the SDW tran- 
sition point, this additional interaction is exponentially 
small compared to a direct pair hopping term. As a re- 
sult, at weak coupling, the pairing does not come from 
spin-fluctuation exchange. Spin-fluctuation mechanism 
may, indeed, become dominant at moderate/strong cou- 
pling, but to verify this one has to extend the present 
approach to uf''^ > 1. 



14 



VII. ACKNOWLEDGMENT 

I acknowledge with many thanks useful discussions 
with Ar. Abanov, E. Abrahams. A. Bernevig, J. Be- 
touras, A. Carrington. P. Coleman. V. Cvetkovic, D. 
Efrcmov, I. Eremin, L.P. Gorkov, K. Haulc, P. Hirsh- 



fcld, Jiangping Hu, H-Y. Kee, Y. B. Kim, M. Korshunov, 
G. Kothar, D-H. Lee, J-X. Li, D. Maslov, I Mazin, J-P. 
Paglione, D. Podolsky, Ph. Phillips, R. Prozorov, D. 
Scalapino, J. Schmalian, Z. Tesanovic, M. Vavilov, A. 
Vishwanath, and A. Vorontsov. This work was supported 
bv lnsf-dmr 06044061 



^ M.R. Norman, Physics 1, 21 (2008); C. Xu and S. Sachdev, 

Nature Physics 4, 898 (2008). 
^ Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. 

Am. Chem. Soc. 130 3296 (2008). 
3 X.H. Chen et al, Nature 453 761 (2008). 

* G.F. Chen et al, Phys. Rev. Lett. 100 247002 (2008). 
^ Z.-A. Ren et aZ.larXiv:080 3.4283l . 

^ M. Rotter, M. Tegel and D. Johrendt, arXivxond- 

mat/0805.4630 (2008)^ 

X.C.Wang et al, arXiv:0806.4688k 3 (unpubished). 

* F.-C. Hsu et ai, Proc. Nat. Acad. Sci. 105 14262 (2008). 
^ J. Dong et al, Europhys. Lett. 83, 27006 (2008). 

^° C. de la Cruz et al. Nature 453, 899 (2008). 
" T. Nomura et al, arXiv:0804.3569 . 

H.-H. Klauss et al, Phys. Rev. Lett. 101, 077005 (2008). 
" Q. Si, and E. Abrahams, Phys. Rev. Lett. 101, 076401 

(2008). 

" M. Daghofer et al, Phys. Rev. Lett. 101, 237004 (2008); 
A. Moreo, M. Daghofer, J. A. Riera, E. Dagotto, 
larXiv:0901.3544l 

F. Ma, and Z.-Y. Lu, Phys. Rev. B 78, 033111 (2008). 
" M. M. Parish, J. Hu, B. A. Bernevig, Phys. Rev. B 78, 

144514 (2008). 

Q. Si, E. Abrahams, J. Dai, J-X Zhu, "arXiv:0901.4112l 
S. Lebegue, Phys. Rev. B 75, 035110 (2007); M.M. Kor- 
shunov, and I. Eremin, Phys. Rev. B 78, 140509(R) (2008). 
D. Singh and M.-H. Du, Phys. Rev. Lett. 100, 237003 
(2008). 

^° L. Boeri, O.V. Dolgov, and A.A. Golubov, Phys. Rev. Lett. 
101, 026403 (2008). 

LL Mazin, D.J. Singh, M.D. Johannes, and M.H. Du, Phys. 
Rev. Lett. 101, 057003 (2008). 
2^ K. Kuroki et al, Phys. Rev. Lett. 101, 087004 (2008). 
K. Haule, J. H. Shim, G. Kotliar, Phys. Rev. Lett. 
100, 226402 (2008); J. H. Shim, K. Haule, G. Kotliar, 
arXiv:0809.0041. 

C. Liu et al ., Phys. Rev. Le tt. 101, 177005 (2008); T. 
Kondo et al, 'arXiv:0807.0 815l . 

2^ D.V. Evtushinsky et al., ar Xiv:0809.4"455l . 

D. Hsieh et al, arXiv:0812.2289l . 
H. Ding et al., arXiv:0812. 05341 . 

2* A.L Coldea, et.al., Phys. Rev. Lett. 101, 216402 (2008). 

S. E. Sebastian et al, J. Phys.: Condens. Matter 20 422203 

(2008). 

^° J. Yang et al, arXlv:0807.1040 

A. Chubukov, D. Efremov and I. Eremin, Phys. Rev. B , 

78, 134512 (2008). 

M.T. Rice, Physical Review B 2, 3619 (1970). 

N. Kulikov and V.V.Tugushev, Sov. Phys. Usp. 27, 954 

(1984), [Usp.Fiz.Nauk 144 643 (1984)] and references 

therein. 

L.V. Keldysh and Yu.V. Kopaev, Sov. Phys.-Sol. State 
Phys. 6, 2219 (1965); J. De Cloizeaux, J. Phys. Chem. Sol. 



M.G. Vavilov, and A.V Chubukov, 



26, 259 (1965); B.L Halperin and M.T. Rice, Sol. State 
Phys. 21, 125 (1968). 

V. Cvetkovic and Z. Tesanovic, arXiv.org:0808.3742 
(2008). See also V. Stanev, J. Kang, and Z. Tesanovic, 
Phys. Rev. B 78, 184509 (2008). 

V. Barzykin and L.P. Gorkov, JETP Lett. 88, 142 (2008). 
S. Graser, T. A. Maier, P. J. Hirschfeld, D. J. Scalapino, 
l arXiv:0812. 0343l. 

S.-L. Yu, J. Kang, and J.-X. Li, arXiv:0901.0821 . 
M.M. Korshunov and I. Eremin, Phys. Rev. B 78, 
140509(R) (2008) 
A.B. Vorontsov, 
larXiv:0812.2469] 

L. Boeri, O. V. Dolgov, A. A. Golubov, Phys. Rev. Lett. 
101, 026403 (2008). 

Y.Yanagi,Y. Yamakawa, and Y. Ono, J. Phys. Soc. Jpn. 
77 123701 (2008); H. Ikeda, J. Phys. Soc. Jpn. 77, 123707 

(2008) 

Y. Bang, H.-Y. Cho i, and H. Won,[ar 2Civ: 0808.3473l . 

Fa Wang et al, larXiv:0807."0498l F. Wang et al, 

arXiv:0805.3343 ; F. Wang et al, arXiv:0807.0498 

K. Seo, B. A. Bernevig, J. Hu, Phys. Rev. Lett. 101, 206404 

(2008) 

In the diagrammatics, RPA corresponds to summing up 
separately ladder series of diagrams and bubbler series of 
diagrams. To study the flow of the couplings one has to 
add renormalizations in the direction transverse to either of 
these channels, i.e., perform parquet-type RG calculations 
(see text). 

T. Y. Chen et al. Nature 453, 1224 (2008). 

A.D. Christians on et al, N ature 456, 930 (2008); M.D. 

Lumsden et aZ.. larXiv:0811. 4755 . 

T.A. Maier and D.J. Scalapino, Phys. Rev. B 78, 
020514(R) (2008). 

Y. Nakai et al, J. Phys. Soc. Jpn. 77, 073701 (2008); Y. 
Nakai et al, 'arXiv:0810.3569 . 

K. Matano et al, Europhys. Lett. 83, 57001 (2008). 
H.-J. Grafe et al., Phys. Rev. Le tt. 10 1, 047003 (2008). 
R.T. Gordon et al, arXiv:0810 T2295l ; R.T. Gordon et 
al., arXiv:081 2.3683l : J.D. Fletcher et a/., larXiv:0812.3858l 
(unpbulished). 

K. Hashimoto et al., Phys. Rev. Lett. 102, 017002 (2009). 

D. Parker et al, Phys. Rev. B 78, 134524 (2008). 

Y. Senga and H. Kontani, J. Phys. Soc. Jpn. 77, 113710 

(2008) 

A.B. Vorontsov, M.G. Vavilov, and A.V Chubukov, 

larXiv:0901.07l9l 

P. Ghaemi, F. Wang, and A. Vishwanath. larXiv:0812.0015l 
(unpublished) . 

D. Parker and L Mazin, arXiv:0812.4416' (unpublished); J. 
Wu and Ph. Phillips, arXiv:0901.0038 (unpublished). 
A.V. Chubukov, I. Eremin and M.M. Korshunov, 
.arXiv:0901.2102. 



15 



Cenke Xu, Markus Mueller, and Subir Sachdev, Phys. Rev. 
B 78, 02 0501(R) (2008); Cenke Xu, Yang Qi, and Subir 
Sachdev, larXiv: 0807. 15421 

P. Chandra, P. Coleman, and A.I. Larkin, Phys. Rev. Lett. 
64, 88 (1990). 

see. e.g., C. Cao, P.J. Hirschfeld, H.-P. Cheng, Phys. Rev. 
B 77, 220506(R) (2008). 

Note that our two-band model differs conceptually from 
the two-band model based on the dyz and dxz orbitals 
considered by S. Raghu et al [Phys. Rev. B 77 220503(R) 
(2008)]. Band structure calculations^li show that not only 
dyz and dxz orbitals, but all Fe 5d-orbitals contribute to the 
bands crossing the Fermi level, and, besides, Fe 5d-orbitals 
strongly hybridize with the As p-states. The tight-binding 
model of the pnictides then should include 8 bands in the 
unfolded BZ. Our model with two hybridized orbitals and 
two bands crossing the Fermi level is not related to actual 
Fe 5d and As p -orbitals, and is simply a phenomenological 
minimal model to describe SDW and superconductivity. 



D. Podolsky, H-Y. Kee, Y. B. Kim, larXiv:0812 . 2907 
A.T. Zheleznyak, V.M. Yakovenko, and I.E. Dzyaloshin- 
skii, Phys. Rev. B 55, 3200 (1997). 

W.L. McMillan, Phys. Rev. 167, 331 (1968). 

J. Wu, P. Phillips, and A.H. CastrcHNeto. larXiv:0805. 21671 

P.Fulde and R.Ferrell, Phys. Rev. 135 A550 (1964); A. I. 

Larkin and Y.N. Ovchinnikov, Zh. Eskp. Teor. Fiz. 47 1136 

(1964) [Sov. Phys. JETP 20, 762 (1965)]. 

A. M. Tsvelik, private communication 

see e.g., M. Eschrig, Advances in Physics, 55, 47, (2006) 

and references therein. 

A.V. Chubukov, D.L. Maslov, S. Gangadharaiah, and L.I. 
Glazman, Phys. Rev. Lett. 95, 026402 (2005); Phys. Rev. 
B 71, 205112 (2005); A.V. Chubukov and D.L. Maslov, 
unpublished. 

see, e.g., N. Shah and A. Lopatin, Phys. Rev. B 76, 094511 
(2007) and references terein. 



